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I. INTRODUCTION 

The depinning transition of a slowly driven elastic in- 
terface in systems with quenched disorder is recognized 
as one of the paradigms of non-equilibrium critical phe- 
nomena. As the driving field is slowly increased and fine- 
tuned to a critical value, the elastic interface undergoes a 
transition from a pinned state to a moving regime. Near 
this pinning-depinning transition, the motion of the inter- 
face proceeds in spatio-temporal avalanches with scale- 
invariant statistics. Avalanches are observed in a variety 
of dynamical systems, including the Barkhausen noise in 
soft magnetic materials [3 the charge density wave 
depinning [3l |4] , the motion of vortices in superconduc- 
tors [5^ , the seismic activity in earthquakes [6] , the acous- 
tic emission in mesoscopic crystal plasticity [7] , and frac- 
ture propagation [5]. Although these systems are very 
different in microscopic detail, the large-scale statistical 
properties of their avalanche fluctuations appear to be 
universally captured by the mean field theory of the de- 
pinning transition [U [HJ [5HT^ . 

An avalanche is a burst of forward motion described 
by a velocity profile v(x,t), where x labels a point along 
the interface. The motion begins from rest [v — 0) at 
time t = Q and ends when the system comes to rest 
again for the first time, at time T. An elastic interface 
driven at a slow, constant rate through a disordered 
medium proceeds forward in a series of avalanches with 
durations that are distributed according to a power law 
P{T) ^ j^-a jj-^ -(-j^j, limit of weak elastic coupling and 
infinite system size Another power-law distributed 
quantity of interest is the global avalanche size, given by 
S — J v{x, t)d'^xdt, which is the total forward mo- 

tion of the interface center of mass, or, in the domain 
wall picture, the total change in magnetization during 
the avalanche [U [TOl US]- Here, we will focus on an- 
other quantity: the maximum collective velocity during 
an avalanche Vm = maxo<t<T J v{x,t)cfix. 

Several acoustic emission (AE) experiments on small- 
scale plastic deformations of single crystals, e.g. ice, Cd, 



Zn, and Cu, report robust power law scaling in the statis- 
tics of the maximum amplitude A,„ of acoustic waves 
emitted during plastic slip avalanches. The probability 
distribution density of the maximum AE amplitude fol- 
lows a power-law tail P{Am) ~ ^m'^i with an exponent 
/i ~ 2 d4.-18j . Since many slip avalanches are required 
to obtain good statistics for P{A^), the deviations in the 
values of /i could depend on the experimental resolution. 
Nevertheless, the exponent is remarkably robust to varia- 
tions such as loading mode, type of crystal, temperature, 
forest hardening effect, or plastic anisotropy [TBHIH]- Un- 
der certain conditions, it is argued that the maximum 
amplitude Am is a measure of the area swept by the fast- 
moving dislocations |14[ I15j . The same power law ex- 
ponent for the maximum velocity distribution has been 
observed using high- resolution extensometry |17) . 

The picture that emerges from the AE experiments is 
that crystal deformations due to the intermittent motion 
of dislocations are characterized by seismic events analo- 
gous to those in the dynamics of plate tectonics. Crystal 
plasticity as a critical phenomenon has also been con- 
firmed by other experimental techniques such as high- 
resolution extensometry |17| and compression of nano- 
pillars and micro-pillars |19[ I20j , as well as by numer- 
ical simulations of discrete dislocation dynamics [7| [^Tl - 
[53]. Moreover, the statistics of slip avalanches-i.e. the 
avalanche size and duration distributions as well as the 
power spectrum-seem to agree remarkably well with the 
mean field theory of interface depinning 23-26 . 

Recently, in Ref. [21] we derived the value of the expo- 
nent n from a mean field theory of the depinning tran- 
sition. In this paper, we present the calculations that 
led to this result in more detail and extend the results 
to include the effect of a nonzero driving rate and cutoff 
scaling functions. Although motivated by experiments on 
crystal plasticity, our results are expected to be broadly 
applicable to other systems that fall in the mean field in- 
terface depinning universality class, such as earthquakes 
and Barkhausen noise in soft magnets. In Table |l] we 
give a summary of the predictions of the mean field the- 
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Quantitiy Form MFT values 

Avalanche size distribution P{S) ~ S"^ Tik^'" S) r = |, cr = | 

Duration distribution P{T) ~ T~"g\k''''T) a = 2, vz = l 

Maximum velocity distribution P{vm) ~ Vm^'H{k''vm) M = 2, p = 1 
Maximum energy {E^ = ^i) distribution ~ E~^^''^^^''^'''H{kP eU'^) f + | = |, 2p = 2 

Stress- integrated size distribution -Pint(>S') ~ t + a = 2 

Stress- integrated duration distribution -Pint(T) ~ T^''^^^' q -f ^ = 3 

Stress-integrated maximum velocity distribution Pint('Um) ~ Vm '^ '' '' /i + ^ = 3 

Stress-integrated max energy (E^ = v'm) dist. Pint{Em) ^ ErS^'^^ '''' ^ ^ l(/^+p) + l~2 

Average size vs. duration {S\T) ~ avz = | 

Average max velocity vs. duration {vm\T) ~ ^ = ~ 1 

Average max energy {Em = vs. duration {Em\T) 2^ = 2^ = 2 

Average max velocity vs. size {vm\S) ~ 5"^'' ap = = | 

Average max energy {Em = Um) vs. size (i5m|<S') ~ g^o-p 2(jp — 2^J^ = 1 

TABLE I. Summary table of the mean-field theory (MFT) exponents and scaling relationships for avalanches statistics in slowly- 
driven (the driving rate c — >■ 0) interfaces near depinning. The first four lines give the power law and cutoff exponents for for 
the size S, duration T, maximum velocity Vm and maximum energy Em = observables. The script letters denote universal 



scaling functions. The first three distributions are shown exactly in Eqs. (731, (72 I, and (78 1 and the fourth follows from a 
simple change of variables. The parameter k gives the distance to criticality for the system; in Barkhausen noise experiments, 
it represents the demagnetizing field while in steady-state plasticity scenarios (including earthquakes) it is proportional to the 
stiffness of the coupling between the system and the driving. In the stress-controlled situation, k can be replaced by Fc — F, 
where F is the external stress and Fc is the critical stress, and the exponent predictions will remain the same. The second 
four lines give stress-integrated exponents and the exponents are the expected outcome of plasticity experiments where the 
external force or stress is increased gradually until failure and avalanches occur along the way. The last five lines give how the 
size, maximum velocity and maximum energy of an avalanche scale with its duration, and then how the maximum velocity and 
energy scale with the size. The first is well-known [1] 1131 1271 128| and the last four can be obtained by taking the averages of 
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ory, including our new results for the maximum velocity. 
The predictions we quote are catered to both Barkhausen 
noise and plasticity experiments (but see the caveats at 
the end of Section VI ) , the latter in both steady-state and 
stress-controlled situations (see the caption to Table |l| . 

Some of the distributions that we derive in this paper 
are relevant to the field of extreme value statistics (EVS), 
which is of importance to many statistical systems, from 
natural disasters, e.g. floods [30], earthquakes [31], snow 
avalanches [35], to financial crises [331 ^^'^ social reac- 
tions to extreme events [33]. Most of these natural sys- 
tems have a complex evolution characterized by catas- 
trophic, extreme events that appear at intermittent inter- 
vals and have a large impact on the long-term behavior 
of the systems |35j . In physics, determining the likeli- 
hood of these rare events is important to a variety of 
problems. For example, the low-temperature thermody- 
namics of spin glasses is governed by the statistics of 
low-energy states [SG^; the velocity of traveling fronts is 
selected by the extreme value for a wide class of ini- 
tial conditions similar to random bisections that occur 
during fragmentation processes |37]; the level density of 
a noninteracting Bose gas has an asymptotic behavior 
governed by the three limit laws of classical EVS [38] 
due to a deeper connection between the sum of corre- 
lated random variables and extreme values [39] : much 
recent progress on EVS in systems with strong correla- 



tion has been achieved in the context of Gaussian height 
fluctuations of kinetically growing interfaces [3DH33) ; ex- 
treme statistics applies to the largest/smallest eigenvalue 
of random matrices, also corresponding to the extreme 
positions occupied in a Coulomb gas with logarithmic 
interactions [35] . 

Classical extreme value statistics applies to indepen- 
dent, identically distributed random variables {xi}fLi 
with the parent distribution p{x). Under these assump- 
tions, the distribution that the extremum (here taken 
as the maximum) is less than a given value to, namely 
C{m) = Prob(maxi<i<Ar < to), approaches in the 
limit of iV — oo one of the three types of asymptotic 
laws depending on the tail of p{x): i) when p{x) decays 
faster than any power law and is unbounded, e.g. a Gaus- 
sian or an exponential function, the EVS is governed by 
the Fisher-Tippett distribution C(m) ^ exp(exp(— to)); 
ii) when p{x) decays as a power law p{x) ~ x~°', with 
a > 0, the extreme value distribution approaches the 
Frechet limit function, C{m) ~ exp(— to~"); iii) the 
WeibuU distribution function C{m) ~ exp(— (— to)"), 
with a > 0, determines the EVS class for the random 
variables that are distributed on a bounded interval. Al- 
though a general theory of extreme value statistics for 
arbitrary correlations and parent distribution p{x) has 
yet to emerge, substantial progress has been made in 
several cases where the EVS is exactly solvable for cor- 
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related variables. For instance, Berman's theorem says 
that the EVS of a weakly-correlated Gaussian process is 
also governed by the Fisher-Tippett distribution [JSJ ST] . 
This corresponds to a power spectrum density that de- 
cays with the frequency as Sf ^ /~°, with an expo- 
nent < a < 1. Indeed, the EVS of time records 
with long-term persistence of Gaussian distributed fluc- 
tuations with a < 1 converges to the Fisher-Tippett dis- 
tribution [H]. The same asymptotic law also determines 
the distribution of maximum heights of periodic, Gaus- 
sian l//-noise |431 149| . where the maximum is measured 
relative to the mean. 

On the other hand, the extreme value statistics of a 
correlated Gaussian process with a > 1 typically has 
a simple scaling form with the duration, with a scaling 
function strongly dependent on several parameters, such 
as the boundary conditions, the value from which the ex- 
tremum is measured, as well as other ordering constraints 
on the time evolution [42l |43]. For example, different 
scaling functions are obtained for the maximum heights 
of periodic Gaussian interfaces: if the maximum is mea- 
sured relative to the spatially averaged height, the cor- 
responding EVS is determined by the so-called Airy dis- 
tribution function |41H44| . whereas measuring the maxi- 
mum relative to the boundary value leads to the Rayleigh 
distribution [501 EI] • We will discuss the extreme value 
distribution for an avalanche in a mean field theory of 
interface depinning known as the Alessandro-Beatrice- 
Bertotti-Montorsi (ABBM) model [52] and will show 
that this avalanche signal can be viewed as a sequence 
of strongly-correlated, non-identically distributed, non- 
Gaussian variables. The probability distribution func- 
tion (PDF) of the maximum velocity inside avalanches 
of fixed duration T follows a universal scaling form 
P{v^\T) = {2vmTy'^/^F{y/2v^/T), with a scaling func- 
tion F(x) that can be derived exactly by a mapping to 
an equivalent problem of random excursions of Brownian 
motion in a logarithmic potential. 

This paper is structured as follows: In Section II, we 
review the ABBM model of interface depinning and de- 
rive it from a discrete mean field theory. Then in Section 
III, we outline our methods and use them to demonstrate 
a derivation of the velocity distribution at a given time 
in an avalanche of a given duration. We also comment 
on the relationship between the Ito and Stratonavich in- 
terpretations of the square-root multiplicative noise in 
the model, in particular why they both predict the same 
parabolic shape for the average velocity as a function of 
time inside an avalanche. Following that, in Sections IV 
and V, we use methods similar to those of Section III to 
derive the distributions of peak velocities for avalanches 
of given durations and sizes. The main results for fixed 
durations at the critical point have been recently re- 
ported in Ref. [53]. Here, we present the details of the 
calculations and extend the work to also include the ef- 
fects of nonzero sweep-rate c and the maximum velocity 
statistics for fixed avalanche size. Finally, in Section VI, 
we derive the peak velocity distribution integrated over 



all avalanches, as well exact scaling functions describing 
the tails of size and duration distributions. Concluding 
remarks and a summary are provided in Section VII. 



II. ABBM MODEL 

Mean-field avalanches in the motion of interfaces near 
depinning have been extensively studied in the realm of 
the ABBM model, which was proposed by Alessandro et 
al. inRefs. |52j- The model was initially derived to mimic 
the dynamics of a single domain wall in soft magnets ex- 
hibiting Barkhausen noise, but has become a 'standard 
model' for interface depinning in the presence of long- 
ranged forces pQ [ini HH [S3] . The dynamics of a single 
interface propagating across a disordered material results 
from a competition between internal elasticity of the in- 
terface, interactions with the quenched disorder experi- 
enced by the moving front and an external driving field 
that is increased at a constant rate. In the mean-field ap- 
proximation, the motion of the center-of-mass position of 
the interface u{t) = J u{x,t)d'^x is given as 

|..(^,-.)+n.o, (1) 

where c/k is the constant pulling rate, k is an elastic 
coupling to the driving force, and the effective pinning 
force F{u) is assumed to be Brownian correlated in u, 
i.e. {F{u)) = and 

{\Fiu)-Fiu')\^)=2D\u~u'\. (2) 

It is important to note that this does not imply long- 
ranged correlations of disorder in the medium. In fact, 
under a physically sensible short-ranged disorder corre- 
lator {F{x,u)F{x' ,u')) ~ 5{x — x')A{u — u'), where A{u) 
is peaked at the origin, the effective correlations in the 
total pinning force F{{u{x)}) = / F{x,u{x))d'^x in the 
center of mass coordinate u are long-ranged because dif- 
ferent pieces of the interface move forward at different 
times and the smaller the fraction that moves, the less 
the total pinning force changes . 

The evolution of the average velocity of the interface, 
defined as v — du/dt, can be obtained by differentiating 
Eq. ([!]) with respect to u. Thus, 

dv , c , . 

-r = -k + - + w{u), (3) 

du V 

where v{u) = v{t) and = dF/du, so that 

{wiu)w{u')) = 2D6{u - u'). (4) 

Eq. ([3| can be interpreted as the motion v{u) of a 
Brownian particle in a logarithmic potential U (u) — 
ku — clog(M). We will assume that the motion obeys 
v{x, i) > at all times, so that u(t) is monotonic and we 
can choose to study the system's evolution in either u or 
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t. By multiplying with du/dt on both sides in Eq. (|3]), 
the evolution in time can then be written as 



dv 
dt 



= —ki 



(5) 



where the multiplicative noise r](t) = 
has zero mean, obeys 



'vw is Gaussian, 



(6) 



and is interpreted in the Ito sense. The interpretation 
can be verified by demanding that Eq. ([5| predict the 
correct steady state distribution 52J 



P{v, t — >■ oo) cx V 



-l+c/D 



cxp{—kv/D), 



(7) 



which can also be derived from Eq. ( [3|) fT]. 

Incidentally, we notice that Eq. ([5|) is identical to 
the equation satisfied by interest rates in the Cox- 
IngersoU-Ross bond pricing model [51]. A similar y/v- 
multiplicative noise process also appears in reaction- 
diffusion systems driven by internal noise [551 . The 
dynamics of this process is characterized by power-law 
statistics of avalanche sizes and durations, and long-range 
temporal correlations with power spectrum S{f) ~ /~", 
where / is the frequency and a > 1. 

An avalanche of duration T corresponds to an excur- 
sion of v{t), which is a path that starts and ends with 
v{0) = v{T) = with v{t) > for < t < T. We study 
the statistics of avalanche durations and velocities con- 
ditioned on durations from Eq. ([5|. Alternatively, the 
avalanche size distributions as well as the velocity statis- 
tics conditioned on avalanche sizes can be studied from 
Eq. ([^. In this parametrization, an avalanche of size 

S = Jf^ v{t)dt corresponds to an excursion of v(u) along 
the coordinate u (integrated velocity). This corresponds 
to a path that starts and ends at v{0) — v{S) — with 
v{u) > for < u < S*. 



A. Other formulations of the mean field theory 

This mean field theory has two other equivalent con- 
structions that provide useful physical pictures. The first 
uses integer time t — jSt and velocity variables and is a 
Markov process where the number of cells rit+i that fall 
at time step t + 1 is drawn from a Poisson distribution 



-p{nj) 



Ij + l'. 



(8) 



where pijit) — (n^+i). This rule can be derived from 
the following picture, (see e.g. Refs. [51 [5B]): We first 
divide the interface into N cells of equal volume that 
are coupled to one another with a mean-field coupling 
(infinite range) and also each coupled to a driving force 
that pulls the interface at an average velocity Va- This 



means that the total local stress on the i'th cell is given 

by 



n = {J/N)Y^{uj - U,) + K{Vat - 



Ui) 



(9) 



where J is the elastic coupling between the cells, K is the 
coupling to the driving point and the Ui are the positions 
of the cells along the direction of motion. 

The dynamics is that a cell slips forward by an amount 
5ui as soon as its local stress is above a threshold r^^^. 
This raises the stress on every other cell by an amount 
Stj = { J/N)5ui and decreases cell i's stress by an amount 
\6Ti\ = (J + K)5ui. Assuming that the stress drops, 
5Ti, and failure stresses, j, are distributed very nar- 
rowly around the values 5t and r^, respectively, each 
cell's stress will be approximately confined to the interval 
[ts — St, Ts] ■ Then, each cell that did not slip at a given 
time step will experience a stress increase of 



ATf = rit 



J6t 

N{J + K) 



Kvjt 



(10) 



where nt is the number of cells that slipped at time t and 
dt is the duration of a time step. If we further assume 
that the stresses of the cells that did not slip at time t 
are uniformly distributed within [r^ j — Sr^Tg i] , the cells 
that fall at time t + 6t will be the ones with stresses 
greater than r^^i — Art. This means that each cell has a 
Art /St chance of slipping at the next time step t + St 
and that the average number of cells that slip is 



H+dt/ 



N 



At* ^ 

St (J + K) 



NKvJt 



(11) 



Writing the average number of cells that slip per time 
step in terms of the pulling velocity as — VaNSt/Su 
and rearranging gives 



p(nt) = (nt+st) 



Jilt + KUa 

J + K 



(12) 



which has (rit+st) — {i^t) = Ua as a, steady state aver- 
age. Since we expect an average of p{nt) cells to slip in 
the next time period, rit+st is Poisson distributed as in 
Eq. ([8|. 

As also pointed out in Ref. [57, an update rule like 
Eq.(|8]) derived in the shell model of the random field 
Ising model [57] is equivalent to the ABBM model in 
the continuum limit. This follows from the Gaussian ap- 
proximation to the Poisson update rule, which is valid 
for large p{nt) and is given by 



nt+st = pint) + ^fp(nt)m 



(13) 



where rjt is a univariate Gaussian random variable with 
zero mean. Assuming that K/ J is small and that terms 
multiplied by it can be neglected in the fluctuation term 
when /o(rif) is large, we have 



nt+6t - nt 



K 

7 



{fit - Ha) + ^/THr|t 



(14) 
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The number of cells that slip during a time step is re- 
lated to the instantaneous center of mass velocity by 
Vt — ntSu/{NSt), so, multiplying through by the con- 
version factor and taking the continuum limit in time 
St — >■ 0, we arrive at the ABBM equation 



v{t) 



v{t)) + ^2Dv{t)i{t). 



(15) 



K 



c = kva, rjit) = rit/y'6t, and 2D 



where k = -jg^, u, =. n,ua, — in/ v ui,, cui^ = j^j^ 
The form of the discrete time equation shows that the 
Ito interpretation is correct. Therefore, this model is 
equivalent to the ABBM model at large velocities. Notice 
also that for K — 0, the continuum limit of the discrete 
model's steady state equation 



P(n) = ~ / dn' 



,(n')"e-" 



P(n') (16) 



is satisfied by the ansatz P(n) = n ^, which agrees with 
the c,k^Q limit of Eq. 

The ABBM model is also the continuum limit of a 
point process in which the dynamics proceeds in discrete 
jumps with random waiting times Tj — tj — tj-i in be- 
tween that are power-law distributed (see e.g. Ref. [SS] 
and references therein). In particular, the discrete pro- 
cess that leads to Eq. Islat the critical point has f (r,) ~ 

For a stationary Poisson process with a constant mean 
A, the waiting time distribution is P{t) — Ae"'^'^, 
whereas for a modulated Poisson process with A a 
stochastic variable with a probability distribution -P(A), 
it is 



P(r) 



(17) 



The steady state probability of the jumps is P{n) ^ n^^, 
so P(A) ~ A~^, which in turn implies that P{t) ^ r^^ 
in the continuum. This is consistent with the rule 
P{Tj) ~ in the discrete case since the probability 
of a random time being in a waiting interval of duration 
T is proportional to t. Since Eq. ^ can be derived as the 
continuum limit of a sequence of discrete pulses, the Ito 
interpretation is again seen to be appropriate. 



B. Numerics 

In numerical simulations, we use an Euler - Maruyama 
method to integrate Eq. ([5| with D = 1/2. This scheme, 
however, does not preserve the positivity constraint of 
the solution. Nonetheless, since we are interested in 
avalanches, i.e. the evolution in between zero crossings, 
we discard the times when the velocities become neg- 
ative. An avalanche is started from v = 1 and evolved 
until u < at which time it is declared over and recorded. 
Then a new avalanche is restarted from v = 1 with a dif- 
ferent realization of the noise. Since the noise is uncorre- 
lated, this closely resembles the reflecting boundary con- 
dition at D = in the steady state evolution of the model. 



This method is simple and fast, but the disadvantage is 
that the crude integration rule and boundary conditions 
introduce inaccuracies for small avalanches (T <C lO'^, 
S ^ 10^). Another method tried was directly imple- 
menting Eq. Q which obtained results consistent with 
the first method for k = c = 0, with a somewhat faster 
convergence to the continuum limit seen in the scaling 
functions. 



III. VELOCITY STATISTICS FOR 
AVALANCHES OF FIXED DURATION 

With the change of variable x = 2y/v and Ito's lemma, 
Eq. ([5]) can be transformed into an additive noise equa- 
tion 



dx 



k 

2" 



2c -D 



(18) 



The problem is thus mapped to an over-damped Brown- 
ian motion of a particle confined to the right half plane 
in a potential V{x) = ^kx'^ — (2c — _D)log(a;). Near 
criticality, fc — >■ and the problem reduces to Brown- 
ian motion in a logarithmic trap. For an initial condition 
x(0) — 0, the particle may execute an excursion with 
x{t) > for a duration T until its position returns to the 
origin for the first time. This excursion corresponds to 
an avalanche in the v{t) during the interface propagation. 
In the following, we calculate the distribution of the in- 
stantaneous displacement during an excursion of size T, 
namely P(x,t\T), and, by a corresponding transforma- 
tion of variables, the conditional distribution of velocities 
in avalanches of fixed durations, P{v,t\T). This allows 
us to determine the avalanche shape. 

Following a standard path-integral method, the proba- 
bility distribution for the process defined by Eq. ( 18 ) can 
be written as a path integral 



P{x,t)<x J Vv{t)5{x - y{t)) 



(19) 



where V{x) is as above and the noise average is per- 
formed with respect to the Gaussian distribution P^ cx 
exp(— / £^(tydt) that produces, e.g. Eq. (jej). Follow- 
ing Ref. we demand that the path is an excursion 
of duration T by fixing the endpoints x(0) — x{T) — e, 
where e is some small value, later taken to zero. Also, 
by applying the positivity constraint that x{t) > for 
< t < T, we have that 



P{x,t\T) = lim 



1 



Z,{T) ^(0)^ 
x5(,v{t)-x) n 



Vy{T)e- ^0 <i^L'^(v,y) 

(20) 



0<r<T 



The noise average has been performed, resulting in the 
Lagrangian (2;, ^) = ^ + £?XM)' and the 
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product of theta functions enforces the positivity con- 
straint. The normahzation factor integrates over all ex- 
cursions regardless if they go through x at time t and is 
therefore given by 



px{T)=e 

Z,{T)= Vx{T)e- drLE(x,x) 

Jx{0)=e 



n ©(^w)- 



(21) 



0<T<T 



The equivalent real-time Lagrangian, obtained by trans- 
forming the imaginary time coordinate t to r = —it and 
factoring out a —1 is L{x,x) = —jj^{—ix + V'{x))'^, 
where i = 37 and V'{x) = The canonical mo- 

mentum corresponding to L{x,x) is then P — jg; = 
jjj{—ix + V'{x)), and the real time Hamiltonian H = 
pi — L — Dp^ — ipV'{x). Therefore, the evolution of the 
Brownian particle in a potential V{x) can be determined 
from its quantum analog that satisfies idrip = Hip. Re- 
placing p = —idx and r = —it, the probability t) 
that the Brownian particle is at position x at time t sat- 
isfies the Fokker-Planck equation 



dP d d'^P 



(22) 



which must be solved with an absorbing boundary con- 
dition at the origin P{0, i) = in order to enforce the 
positivity constraint. For fc = 0, Eq. (22) reduces to 



dP 
~dt 



d_ 
dx 



D-2c 



-P 



D 



(fP_ 
dx"^ ' 



(23) 



which can be solved exactly. 

By a generic separation of variables P{x,t) 



e f{x), Eq. (23) is reduced to the eigenvalue problem 



dx"^ 



l-2cdf 
X dx 



E 
D 



1 ~ 25 



/ = 0, 



(24) 



where c = c/D. This is a modified version of the Bessel 
equation and the particular solution that is well-behaved 
at the origin is 



f{x) = x'j, 



l-c 




(25) 



where Ji_g(a;) is the Bessel function of the first kind 
We have assumed < c < 1 since this is the range in 
which the model exhibits avalanche behavior (otherwise 
the origin is inaccessible in the steady state, as can be 
seen in Eq. ([t])). We can now proceed to calculate the 
P{x,t). For this, we interpret the above path integrals 
as a matrix element between two functions defined on 
[0, cx)) given by 

POO 

{h\e-"'^'-*°^\g)^ = h{x)e-"'^*''°'>g{x)dx. (26) 
Jo 



where the Hamiltonian is Hi = —Dd^ — dxVl{x) with 
the potential Vi{x) = — {2c — D)log{x) and an absorbing 
boundary condition at x = 0. We write the eigenfunc- 
tions fk{x) = x'^Ji-c{kx) with the eigenvalues defined 
as Ek = Dk^. Using the inverse Hankel transform of 
g{x)x~'^ — /p fc(ifcJi_5(fcx)gi_g(A;), the above expression 
becomes equivalent to 

/"OO 

(/i|e-^i |.g)„ = / dxh{x) X 
Jo 

1*00 

X / fcdfce"-°'='(*-*°)xVi_£(/ca;)5i_£(fc), (27) 



where 5i_g(fc) = /q°° Ji_c{kx)g{x)x^~'^dx. 

Hence, the probability of being at position x and time 
t during an excursion of duration T, P{x,t\T), can be 
defined using the matrix elements as 



P{x,t\T) = lim 



e-i-O 



-, (28) 



where {(\e-^^'^\e)^ = ^e(r). Thus, from Eq. ^ with 
h{x^ = b{x — z) and g(x) — S{x — e), we have 



{z\e 



-Hit, 



£)vj = I dx5{x — z) I kdkx'^ Ji-c{kx)e 
Jo 



-Dk^t 



X / ''dyJi_i,{ky)S{y - e) 

10 

-Dk^t 



e^-^z^ 



kdke 



(29) 

Ji_£(/ce)Ji_g(fcz), 







and, in the limit of e ^ 1, the above expression can be 
expanded to leading order 



{z\e-"^'\e)^ « ^TT^TwTT^ / e-^''''Ji,,{kz)k-'-^dk 



2i-'^r(2 - 5) J, 

e^2(l-c) z(Dt)-^+^ _^ 

— — e 

2 J 2r{2-d) 



(30) 



where T{z) is the gamma function. 

In a similar way, we obtain the other transition ampli- 
tude, to the leading order in e, as 



ez 



l-c 



2i-T(2 - a) J, 

'z\2(i-£) e[D{T-t)]-^+^ 



V2/ 2r(2-5) ^ ' 

The normalization constant Z^{T) can also be deter- 
mined with h{x) — 5{x — e) and g{x) — 5{x — e), and 
is given by 



^2 -Dk^T 



POO 

Z,{T) = e [Ji_£(fce)]'e 
Jo 

I g 2DT 

2DT J 



kdk 



10 
e 



'l-c 



2DT 

(£/2)3-2g 

(OT)2-£r(2-5)^ 



(32) 
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FIG. 1. (Color online) Numerical results for the instanta- 
neous velocity v{t) at time t averaged over many avalanches 
of duration T and rescaled to collapse to the form of Eq. ( |34[ ) 
with n = 1. Results for three different values of the dimen- 
sionless driving rate c = c/D are shown. (In the numerics, 
the driving rate c is varied while the fluctuation strength is 
fixed at D = 1/2.) The parabolic shapes with c-dependent 
height predicted by Eq. ( |35[ ) for n = 1 match well with the 
numerics. 



where Ii-^{x) is the modified Bessel function of the first 
kind and c < 2 |59j . Combining these final expressions 
with Eq. (28), we arrive at 



P{x,t\T) 



(x/2) 



3-2c 



T 



r(2-5) \Dt{T-t) 



2-c 



X exp 



ADt{T ~ t) 



(33) 



Finally, the distribution of velocities at a given time t 
inside an avalanche of duration T follows by a change of 
variable v = x"^ /A and is given as 



From Eq. ( |34[ ), we determine the profile of the ri-th mo- 
ments to be 



obtained from the exact solutions to the untransformed 
fc 7^ Fokker-Planck equation |52j 



dtP ^ Ddl{vP) + d^iikv - c)P), 



(36) 



which were obtained by Feller [60 . The exact propagator 
with an absorbing boundary condition at w = is given 

by 



P{v,t.,vo,0) 



5 exp 



j?("+"oe 
1— cxp( — fci) 



-fet 



1 - e-*^* 

(c-l)/2 



(37) 



h-c 



VqVC 



kt 



D{e'^^ - I) 



assuming that < c < 1. Repeating the steps above gives 

\ 2-c 



P{v,t\T) 




D 



(1- 



-kt 



) (1 - e-'^'(T-*)) 



yv (1 - e-^-^) 



-fet 



)(1 



o-k{T-t) 



(38) 



and 



{v-{t)\T) = 



r(n + 2 - 5) 

r(2 - 5) 



D I 



-kt\ 



fc(l 



o-fcT 



(39) 



which reduce to Eqs. (34 1 and (35) when A: — >■ 



Interestingly, the same avalanche shape is obtained 
when Eq. ([s]) is interpreted in the Stratonovich sense, as 
derived in Ref. (53]. In this interpretation, the x = 2^ 
transformation of the equation of motion is 



dx 
'dt 



2c 

X 



(40) 



which is a free Brownian motion in the limit k — c — 0. 
Using the same method as above, the conditional distri- 
bution of velocities in avalanches of duration T is (as- 
suming k ~ 0), 



-Pstrat.(w,i|r) = 



,1/2-c 



T 



3/2-c 



{v-it)\T) 



Strat. 



with n = 1 corresponding to the avalanche shape {v(t)\T) 
(see Fig. [T]). We notice that the avalanche shape and 
higher moment profiles are independent of the driving 
rate c up to a nonuniversal pref actor. We have also veri- 
fied this numerically as shown in Fig. ([l]). 

This method of obtaining Eq. ( [34| has the advantage 
of being more intuitive thanks to the correspondence 
with excursions and we will find that this method contin- 
ues to be very useful when we compute the statistics of 
avalanche maxima below. However, Eq. (34) can also be rameter shift c 



r(3/2-5) \Dt{T-t) 

and the n-th moments are 

r{n + 3/2-d) (Dt{T-t) 



r(3/2-5) V 



T 



(41) 



(42) 



which is indeed identical to the results from the Ito case 
up to a non- universal, constant pre- factor. In Ref. |53| . it 
is conjectured that the Stratonovich interpretation would 
give the same results as the Ito interpretation for all uni- 
versal quantities, but this is not the case. Notice that 
Eq. ( 40 ) can be transformed into Eq. ( 18 ) by a simple pa- 
c — D /2. Therefore, any quantity that 
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FIG. 2. (Color online) Scaling collapses of numerical re- 
sults of P{vm\T), the distribution of maximum velocities for 
avalanches of duration T, onto the form given by Eq. (54 1 
performed for several values of the dimensionless driving rate 
c = c/D. (In the numerics, the driving rate c is varied while 
the fluctuation strength is fixed at D = 1/2.) The correspond- 
ing scaling functions F{x) given by Eq. ( |55[ l fit very well with 
the numerics. 



has c-dependence will be different in the two interpreta- 
tions. This includes the avalanche size scaling exponent 
T and the duration scaling exponent a which have the 
values r = (3 — c)/2 and a = 2 — c (see Refs. [TJ \2S\ and 
below). The Stratonovich interpretation predictions for 
these exponents are then Tg = 5/4— c/2 and as = 3/2— c. 
The fact that the Stratonovich interpretation gives the 
same parabolic avalanche shapes and higher moment pro- 
files is only a reflection of the fact that those quantities 
do not depend on the driving rate c (except in an overall 
prefactor). 



IV. MAXIMUM VELOCITY STATISTICS FOR 
AVALANCHES OF FIXED DURATION 



every point inside the excursion, it follows that the vari- 
ables x{t) are not identically distributed. Furthermore, 
the individual distributions are non-Gaussian for small 
X, but converge to a Gaussian right tail for large enough 
X. The maximum displacement in an excursion is defined 
as 



AI — max x(t), 

0<t<T 



(44) 



with a probability distribution determined from the prop- 
erties of P{x,t\T). This follows from the fact that the 
probability of the maximum being less than a certain 
value M < Xm is the same as the probability that at ev- 
ery instant during the excursion the displacement x(t) is 
less than Xm- The path integral representing this cumu- 
lative probability is 



CeiXmlT) 



1 



Ze{T) 7x(0) = 



Vx{t)e- ^"^''^^'^'^^ 



\[ Q[x{t))Q{x„,-x{t)), 



(45) 



where Le{x, x) 



0<t<T 



1 

4D 



D-2c\ 



and the product over 



the Heaviside functions selects the excursion paths that 
satisfy the constraint that the position at any time is 
less than Xm- The normalization constant Z^{T), cor- 
responding to the path integral over all excursions, was 
calculated previously in Eq. ( 32 ) . At the end of the cal- 
culation, we take the limit as e — >■ 0, so that the path 
starts and ends at the origin, staying between and Xm 
during an avalanche of duration T. On the other hand, 
Ce{xm\T) is also the probability that M < Xm, namely 



Ce(x„,|r) = / dMP,{M\T), 



(46) 



thus, the PDF P^{xm\T) follows directly by a differenti- 
ation of Ce(xm\T) with respect to Xm- 

Using the quantum analogue, the path integral in 
Eq. ( 45 ) is equal to a transition amplitude between the 
position eigenstate |e) and itself, and thus can be com- 
puted by an expansion in energy eigenfunctions: 



We now determine the extreme value distribution of 
the maximum velocity inside avalanches of fixed dura- 
tion T, for fc = and nonzero c. Using the mapping to 
the random excursions, we first determine the statistics 
of the maximal displacements in excursions over a fixed 
interval. Let us consider an excursion {x{t)}t in the time 
interval Q < t < T. The probability distribution at a 
fixed time t is given by P{x,t\T) from Eq. (33 1, which 
can be written 



P[x,t\T) = -Qre 



X 

2^t 



(43) 



wi th Q REjy) = 



y3-2c^-y ^Y{2 



c) and CTf = 
t)/T. Notice that because at is different for 



(47) 



with 



{e\e-"'' \e)b ^ / dxSix - e)e-"'' S{x ~ e), (48) 
Jo 

where the Hamiltonian Hi has square-well (box) bound- 
ary conditions at a; = and x = Xm- With these bound- 
ary conditions, the Hamiltonian has a discrete spectrum 
of eigenfunctions given by 



fn{,x) = X^Ji- 



\lX 



(49) 
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with eigenvalues En = DX^/x^, where A„ is the n'th 
zero of the Bessel function Ji_g(a;). The (5-function in 
Eq. (|48|) can be expanded in this basis as 



S{x ~ e) 



2ei 



(50) 



1.6 
1.4 
1.2 



with /„(a;)'s defined in Eq. (49 ). Inserting this expansion >^ 
into the definition of the matrix element from Eq. ( 48 1 , ^.c^ q , 
we obtain 



|e)b = 2e 



OO 

E 



Jl-. 



(51) 



CO 
8 0.6 



0.4 



Finally, inserting the solution of (T) from Eq. ( 32 ) into 
Eq. (47), and taking the limit of e — ^ 0, we arrive at the 



expression for the extreme value statistics of the random 
displacements C{Xm\T) = hm^-^o Ce(a;m|T), namely 



0.2 




(a) 



S=10 
S=2x 10*^ 
S=4x lO*" 
- Theory 



0.6 0.1 



C{xn,\T) 



OO 



r(2 - 5)x 



4-2c jL^ 



.2(l-c) 



1.2 

V (2DS) 



1.4 1.6 

1/2 



2.2 



.(52) 



The distribution P(xra\T) of the maximum dis- 
placements follows from the cumulative distribution as 
P{xm\T) = d^,^C{xm\T). The corresponding PDF of the 
maximal avalanche velocities is obtained by the change 
of variable Xm ~ '2,^/Vm and the transformation 



P{Vm\T) = P{Xm\T) ^ 

It follows that the PDF P{vm\T) is given by 

1 



k=10 
A k=5xl0" 



Theory 



P{Vra\T) = 



-.F 



2Vn 



V2Dv„,T y \ DT 
with the scaling function given by F(x) ~ 



E 



a: 



2-2c 



X2 



(4 - 25) e 



r(2-5) 0:5-22 [J2_,(A„)]: 

(55) 

We have numerically verified Eq. ( 54 ) for different val- 




ues of c and fc = 0. Our collapsed distributions as shown 
in Fig. [2] are in excellent agreement with the analytically 
predicted scaling function from Eq. (55 1. 



V. MAXIMUM VELOCITY STATISTICS FOR 
FIXED AVALANCHE SIZE 

Using the same technique as in the previous section, 
we determine the distribution of maximum events in 
avalanches of a given size. Here, we consider Eq. ([3| 
with c = and nonzero k, which corresponds to a biased 
Brownian motion 



FIG. 3. (Color online) (a) Scaling collapse of P{vm.\S), the 
distribution of maximum velocities for avalanches of size S, 
onto the form of Eq. (66 1. The rescaled distributions ap- 
proach the analytical scaling function from Eq. ( |67[ ) for suffi- 
ciently large avalanche sizes, (b) Scaling collapse for the over- 
all maximum velocity distribution, P{vm) for quasistatic driv- 
ing (c = 0) and several different values of k (see the caption 
of Table |l] for the physical definition) to the form of Eq. (78 1. 
The scaling function plotted is F(x) — Ca;'^ sinh~^(a;) where 
the constant C — 1.44 is adjusted to fit near the tail and is 
proportional to that given in Eq. ( |78[ ) with D = 1/2. Adjust- 
ment of this non-universal constant factor is required since 
the part of the distribution near the origin that deviates from 
scaling alters the overall normalization. The inset shows the 
same maximum velocity distributions before rescaling. 



dv 
du 



-k + w{u), 



(56) 
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where (w(m)w(m')) = 2DS{u — u'). 

The cumulative distribution C{vm,\S) conditioned on 
fixed avalanche sizes is defined by the path integral 



C,{Vm\S) 



1 



{i(S)=e 



^iW) Jv{0)=e 

H Q{v{u))e{vm-v{u)), (57) 



0<u<S 



where the normalization factor is the path integral over 
the unconstrained excursions, 



pviS)=e 
Ji}{0)=e 



(58) 



0<n<S 



Instead of duration T as in the previous calculation, the 
value of u at the end of an avalanche is the slip size 
of the avalanche S — v{t)dt. The path integral in 
Eq. ( 58 ) is equivalent to the transition amplitude between 



the position eigenstate |e) and itself, and thus can be 
computed by an expansion in energy eigenfunctions: 



(59) 



where H2 — -~Dd^ — kdy with an absorbing boundary 
condition at w = 0. The eigenfunctions for this Haniilto- 
nian are 

Uv) = ^Jlsmiqv)e-''^/'^, (60) 

with eigenvalues satisfying Eg = D (^(^^)^ + q'^^ for < 
q < 00. Therefore, the Z^{S) from Eq. (591 becomes 



Z,{S) = / did{i-e)e 
Jq 



-H2S 



6iv-e). (61) 



We expand the (5-function in terms of the eigenfunc- 
tions as S{v — e) — dqg[q)e^^'"/'^^ sm{qv) where 

g{q) = ^/l dve'^'"/'^^ 5{v — e)sui{qv). Applying this 
to the second delta function and Taylor expanding in the 
lowest order of e, we arrive at 



Z,{S) 



'dgsin2(q,)e-^((^)'+'')^ 



e iD . 



TT (2DS')3/2 



(62) 



Similarly, the numerator path integral in Eq. (57) is de- 



termined by expanding in the discrete set of eigenfunc- 
tion of the Hamiltonian H2 with square-well boundary 
conditions at u = and v = Vm- The eigenfunctions are 



with the eigenvalues E,, = D + There- 

fore, we obtain that 



X exp —DS 



k 

2b 



, (64) 



which, in the limit of e — > 0, leads to 



C{v^\S) = 



^{2DS) 



3/2 



E 



X exp 



n— 1 



riTT 

Vm 



(65) 



Interestingly, k has canceled completely out of the final 
answer. 



Differentiating with respect to Vm in Eq. (65), we ob- 
tain the PDF for Vm given by the scaling form 

^<»"'|^' = 7m^^ {7m) ■ 

where Fs{x) is the Brownian excursion scaling function 



— 3 e 2a: 



\J2-K ■^22 



Using the Poisson summation formula |61j 

00 , — 00 

n— 00 * 71— 00 

the scahng function can be written 

Fs{x) = ^(32n4x3 - 24n2x)e-2"'-'^ 

n=l 

and we can see it has the asymptotic form 



(67) 



(68) 



(69) 



(70) 



(32a;3 - 24x)e- 



, X — )■ 00, 



which guarantees that all of its moments are finite. In 
Fig. [3| we present the numerically computed distribution 
of maximal velocities for different avalanche sizes which 
agrees very well with the analytical result. The rescaled 
PDF P(vr,i\S) collapses onto a scaling form as predicted 
by Eq. 



VI. OVERALL DURATION, SIZE AND 
MAXIMUM VELOCITY DISTRIBUTIONS 



2 . 2 ( UTTV ^ „-kv/2D 



In{v) ^ \l — sm- I I e 

Urn. \ V' 



(63) 



The overall duration distribution is proportional to the 
probability that an avalanche returns to the origin at 



11 



(a) 



c/D=0k=10 
c/D=.2k=10"-' 
c/D=0k=10~' 
c/D=.2k=10~' 
c/D=0k=10"' 
c/D=.2k=10"' 
- Theory c/D=0 
Theory c/D=.2 




10 
kT 




FIG. 4. (Color online) (a) Numerical scaling collapse of the 
avalanche duration distribution P{T) for different values of k 
(see the caption to Table |l] for the de finit ion) and the driving 
rate c = {v)/k to the form of Eq. (721. The scaling func- 
tions plotted are Fe{x) = C{c)e^x'^~^e'^ — 1)"''^ which are 
proportional to the ones given in Eq. I\72\ . The non- universal 
constants are fit by eye as C(0) — 3.0 and C{.2) = 3.375. The 
inset shows the same distributions before rescaling with the 
curves for c = 0.2 offset for visibility, (b) Collapse of the size 
distributions for different values of k to the form of Eq. (73 1. 
The theoretical curve plotted is F{x) — C exp{—x/2) where 
C = 0.6 which is, up to normalization, what is given in 
Eq. (73 1 with D = 1/2. The inset shows the size distribu- 



tions before rescaling. 



tween. This means that, from Eq. (371, we have 



P(T) = \imN{e)P{e;T,e,0) 



lim N(e) 

£^0 



2k 



1 



2-£ g(£-l)fcTgl-c 



r(2-£) 



,(71) 



where N{e) is a proportionaUty constant that cannot sim- 
ply be set by normalization since P{e, T; e, 0) ~ 7^-2-i-c 
for T — > and the normalization integral diverges at this 
limit. This divergence must occur because we expect that 
the return time to the origin gets shorter as the starting 
value e gets closer to zero. To get rid of it, we need to 
set a cutoff T* ^ e defined as the minimum duration 
of an observable avalanche. The natural long duration 
cutoff is set by 1/k. Assuming that these cutoffs are 
well-separated, namely that kT* ^ 1, we have 



p{T) = /c(i -s)(fcr*) 



l-c„feT 



1 



okT 



(72) 



for T > T* . This has power law behavior P(T) ^ 
T~°'Q{kT) with exponent a = 2 — c and G{x) a cutoff 
function, as can be seen by taking the limit /c — > 0. The 
exponent agrees with the one reported in Refs. [Tl[55]. In 
Fig. we numerically verify Eq. ( 72 1 for different values 
of c and k. 

The avalanche size distribution can be computed simi- 
larly. Starting from Eq. ([3| with fc = as in Rcf. [28 and 
using the results above for a random walk in a logarith- 
mic potential, one can obtain the exponent r = (3 — c)/2. 
For fc ^ and c = 0, we can use Eq. (62) and normal- 
ize with a cutoff S* ~ as we did with the duration 
distribution, obtaining 



/S* 



P{S) = z^^e^ro'^''^ 



2S'3/2 



(73) 



The tail behavior of the overall distribution P{vm) of 
avalanche maxima can then be derived in the adiabatic 



limit by integrating Eq. (67) against the overall size dis- 
tribution Eq. ( 73 ) , 



P{v,r.\S)P{S)dS 

V2DS 



2vl 



xFs{x)e 2(20") dx. (74) 



The asymptotic behavior of the integrand allows us to 
take the upper limit of the integral to infinity provided 
Vm ^ \/2DS* . The infinite sum can then be integrated 
term by term using the formula 



time T, assuming we use absorbing boundary conditions 
to guarantee that it will remain strictly positive in be- 



(75) 

where (x) is the modified Bessel function of the second 
kind. In this way, the correction factor to ii^^-scaling 
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becomes 



(76) 



The infinite series can be summed using the formula 
^^nexp(— ncc) = 4/ sinh^(x/2), giving 



1 ( kvjn \^ 




Normalizing with a small velocity cutoff v*^ ^ \/2DS* , 
the distribution is 



P{Vm) 



\ 2D ) 



< \ sinh(4^^ 



(78) 



2D 



which is valid for Vm ^ We have already seen in 
Ref. [29] that the driving rate dependence of the power 
law P{vm) ~ Vm'^ at criticality is fi = 2 — c. In Fig. [3]d, 
we confirm Eq. ( 78 1 numerically. 



Although we have derived the exact functional form of 
P(vm), it is important to note that the fj, = —2 expo- 
nent follows essentially from two facts. The first is that 
the scaling forms given by Eqs. (54) and (661 have the 



exponents one would assume from dimensional analysis, 
and the second is that the integral leading to Eq. (78) 



converges when its bounds are taken to and oo. Under 
these conditions, the fluctuations in Vm are distributed 
narrowly about the average value {vm\T) ~ T, and so we 
expect the maximum velocity to have the same scaling 
exponent as the duration. However, there is no guarantee 
that these conditions hold in all cases. 

In fact, this simple scaling by dimensional analysis runs 
contrary to what one might normally expect for an ex- 
treme value distribution. For example, if one wanted 
to make a simple argument, one could ignore corre- 
lations and boundary conditions and crudely approxi- 
mate the avalanche as a sequence of iV = T/St inde- 
pendent, exponentially distributed variables with aver- 
age value {v\T) = AT where 6t is the duration of a 
time step and A is a proportionality constant with the 
same dimensions as D. Then, the probability distribution 
function of the maximum value is given by P{vm\T) — 
N_^-v^/AT _ ^-v„jAT^N-\ rpj^g expectation value 
of Vm can be determined by using the binomial ex- 
nansion P(v \T] - (jv-i)!(-i)'° (fc+i)„„/AT 

The integration corresponding to the first moment can 

be performed term by term, and the average maximum 

takes the form iv \T] - iV ^^N-l (jv-i)!(-i)'' (ATf 
taKes tne lorm [Vml-i- ) — at /^k=Q k\(N-i-ky. {k+iy 

This is equivalent to {v„r\T) = AT {^{N + 1) + 7) where 

*(^) ~ '^'"d^''^'' ^'^^ Euler's constant. As T — > 00, 
this expression scales asymptotically as AT\og{T/St). 
The logarithmic enhancement is due to the fact that the 
random variable has T/6t independent tries to achieve 
a value well above its average. Integrating this against 



the duration distribution P{T) ^ T"^ gives a max- 
imum value distribution with leading order behavior 
Pi'^'m) ^ log(^m)'^m^i a slowcr dccay than the exact 
answer P{vm) ~ '^m^- However, temporal correlations 
change the picture considerably, preventing the average 
maximum value from scaling faster than the average, so 
the actual scaling law is (i;,„|T) ^ {v\T) ^ T which cor- 
responds to P{Vm) ~ 

In Table [l] we summarize the quasistatic mean field 
theory results for the scaling functions and exponents, in- 
corporating the new results for the maximum avalanche 
velocity. However, it should be stressed that these cut- 
off exponents have not been observed in all cases in ex- 
periments and simulations with long-ranged forces, even 
when the power law scaling exponents are mean field. 
In Ref. PJ], it is argued that in the case when the cut- 
off is caused by a demagnetizing field coupled globally, 
with dtu = —k J u{t, x)d?x, instead of the local coupling 
dtu = —ku{x,t), the cutoff exponents are different from 
the mean field predictions in a way that follows from 
dimensional analysis and is consistent with Barkhausen 
noise experiments and long-range simulations. In our 
case it would change the maximum velocity cutoff expo- 
nent from p = 1 to pfc = |. The effects cancel out and 
do not alter the exponent products that appear in the 
last three rows of Table |l] In other cases [5j , deviation of 
the cutoff exponents from the mean field values may be 
caused by hardening, an effect not accounted for by our 
considerations. 



VII. DISCUSSION AND CONCLUSIONS 

We have calculated exact scaling functions and expo- 
nents for the maximum avalanche velocity statistics in 
a mean field approximation of interface depinning. The 
distribution of maximum events in avalanches of fixed 
duration has a robust scaling form with a scaling func- 
tion that depends on the driving rate, while the distribu- 
tion for fixed avalanche sizes is independent of the elas- 
tic coupling constant in the adiabatic limit. The statis- 
tics of maximum velocities in arbitrary avalanches follows 
by integrating the maximum value distribution for fixed 
avalanche size or duration against the size or duration 
distribution. We find that the distribution of peak ve- 
locities in a train of non-overlapping avalanches has a 
scaling regime with power-law exponent /i = 2 — c fol- 
lowed by cut-off regime for Vm > D/k. Although the 
mean field theory captures very well the universal statis- 
tical properties of avalanches, i.e., the size, duration, and 
maximum value of plastic slip avalanches, it still remains 
an open problem as to what extent it can also describe 
other statistical quantities, such as velocity fluctuations 
during plastic deformations. For instance, discrete dislo- 
cation simulations report that the individual dislocation 
velocity probability distribution follows a power law with 
exponent ^ —2.5 [7 , whereas Eq. ([t]) from mean field the- 
ory predicts an exponent of — 1 for the collective velocity 
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distribution. 

Based on our calculations, we expect the w"'' power 
law scaling of the maximum velocity distribution to be 
a very robust experimental observable, even in the case 
of poor time resolution where the true maximum might 
be missed. In fact, even if all an experiment could ac- 
complish was to measure a random velocity within each 
avalanche, the inverse squared power law prediction does 
not change. This can be seen by a simple scaling argu- 
ment. Namely, Eq. ( 54 1 predicts that the average value of 



the maximum velocity scales linearly with the avalanche 
duration as {vm\T) ~ T, thus, in the same way as the av- 
erage velocity {v\T) — {S\T)/T ~ T. Since the maximum 
does not outpace the average as the avalanches get larger, 
we expect the same power-law behavior for the distribu- 
tion regardless of whether we sample the maximum or 
a random point in the avalanche signal. This implies 
that even if the maxima are taken from a low-resolution 
time series, one should see an inverse squared power law 
distribution in the limit of slow driving, and this may 
make the maximum velocity easier to work with than 
the duration in low-resolution experiments. The scaling 
{vm\S) ^ S^/"^ should show similar independence from 
how well Vm is able to be measured. 

However, one thing that can affect the exponent sub- 
stantially is the driving rate, which must be very slow in 
order to see fj, = 2. In fact the variability in the acoustic 
emission experiments f7l [T4HI8] on the value of /i might 
be related to the fact that the experiments are performed 
at nonzero driving rate. If that is the case, we find that 



the deviation from the adiabatic driving rate enters in the 
mean field exponents as /i = 2— c, in a way similar to that 
for the driving-rate dependent scaling exponents a and 
T in the power-law distributions of avalanche durations 
and sizes. Experiments where the power law exponent 
is studied as a function of applied shear rate could test 
this. 

Finally, we have shown that exact results for scaling 
functions in the maximum velocity distributions are easy 
to derive, even including nonzero sweep rates and tuning 
away from criticality. For instance, our results for the 
maximum velocity distributions P{vm\T) and P{vm\S) 
given duration and size, given in Eqs. (|54| and (66) re- 



spectively have simple scaling forms, the second of which 
is independent of the distance from criticality k. Though 
it might be difficult to obtain sufficient statistics to check 
the scaling functions of these distributions experimen- 
tally, the predicted scaling of the conditional moments 
{v^\T) ~ T" and {v^\S) ~ 5"/^ could be checked much 
more easily. This could serve as a starting point for new 
precision tests of universality in the Barkhausen effect 
where high resolution time series are obtainable and there 
is already strong evidence of mean- field behavior |55] . 
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